Layer Analysis of the Structure of Water Confined in Vycor Glass 
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A Molecular Dynamics simulation of the microscopic structure of water confined in a silica pore 
is presented. A single cavity in the silica glass has been modeled as to reproduce the main features of 
the pores of real Vycor glass. A layer analysis of the site-site radial distribution functions evidences 
the presence in the pore of two subsets of water molecules with different microscopic structure. 
Molecules which reside in the inner layer, close to the center of the pore, have the same structure as 
bulk water but at a temperature of ?&K higher. On the contrary the structure of the water molecules 
in the outer layer, close to the substrate, is strongly influenced by the water-substrate hydrophilic 
interaction and sensible distortions of the H-bond network and of the orientational correlations 
between neighboring molecules show up. Lowering the hydration has little effect on the structure 
of water in the outer layer. The consequences on experimental determinations of the structural 
properties of water in confinement are discussed. 



I. INTRODUCTION 



The observation that in many real situations water is 
confined (in a rock, in a cell, in a microemulsion, in an 
ionic channel, in interstellar bodies etc.) explains the 
large interest devoted to structural and dynamical prop- 
erties of water in restricted geometries. During the last 
decade indeed several experiments have been performed 
by using different techniques and different confining sub- 
strates, that have evidenced a rich variety of phenom- 
ena, sometimes even more puzzling than those known for 
bulk water .El Moreover recent theoretical studies have ev- 
idenced for confined water a thermodynamics more com- 
plex than that of bulk water .□ 

In particular the comparison of the properties of bulk 
and confined water that emerges from different exper- 
iments is sometimes misleading. As a matter of fact 
the slowing down of the translational and rotational mo- 
tion upon confinement, evidenced by nuclear magnetic 
resonance (NMR)H and quasi elastic neutron scattering 
(QENS),Q apparently contrasts with the lower average 
density deduced for confined water with respect to the 
bulk state by small angle neutron scattering (SANS)lI. 
Moreover neutron diffraction (ND) experiments on pure 
D2O have been interpreted in terms of an increased num- 
ber of hydrogen bonds (HB) upon confinement)!] while 
experiments exploiting the isotopic H/D substitution on 
water reveal adduction of such bonds and a lower tetra- 
hedral order ,0'El in agreement with X-ray diffraction.B On 
the other hand the analysis of the temperature deriva- 
tives of the radial distribution functions (RDF) of con- 
fined and bulk D2O suggests that the structural rear- 



rangement upon lowering the temperature is the same for 
the two systems!] More recently the idea has been put 
forward that lowering the amount of water inside the con- 
fining matrix is. e ji nu ivalent to lowering the temperature 
in bulk water BlilJO Finally, besides and in partial dis- 
agreement with these interpretations of the experimental 
findings, two experimental evidences arc well established: 

• Confined water can be easily supercooled down to 
lower temperatures compared to the bulk. 

• The lower is the hydration level, the lower is the 
ice nucleation temperature. 

The observation that the density -profile of water in 
restricted geometries is not uniform,E3ii3 whatever the 
substrate- water interaction is, suggests that most of the 
misfits between the experimental observations reported 
above are due to averaging the response of different wa- 
ter layers, that cannot be easily isolated in the exper- 
iments, unless an accurate study as a function of the 
hydration level is performed. Molecular dynamics (MD) 
simulations performed in realistic environments can be 
then of great help in better interpreting the experimen- 
tal data. As a matter of fact the syster 
more widely used by the experimentalistsoOQ'tit c 
to confine water, namely porous Vycor glassJ2j has also 
the advantage of being a siaiBle. molecular system, that 
can be computer simulatedOtEl In the attempt of giv- 
ing an harmonic interpretation of at least the structural 
properties of water in restricted geometries, we have per- 
formed a layer analysis of the microscopic structure of 
water in Vycor glass, at different temperatures and hy- 
drations. 
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In Sec. II we briefly describe the simulation method. 
In Sec. Ill we discuss how the calculation of the pair cor- 
relation functions can be performed by taking into ac- 
count the finite volume effects. In Sec. IV we introduce 
the layer analysis and present the main results. The last 
section is devoted to the conclusions. 



II. MOLECULAR DYNAMICS OF CONFINED 
WATER 







The behavior of water confined in the pores of Vycor 
glass is influenced by both the geometrical effects and 
the interaction with the hydrophilic surfaces. In order 
to account for these effects in computer simulation we 
constructed a single pore of silica which reproduces the 
most relevant average features of the pores of real Vycor 
glass. As described in more detail in a previous workli3, 
we carved a cylindrical pore of 40 A diameter in a cu- 
bic cell of silica glass obtained by the usual procedure of 
MD simulation. Then we prepared the cavity surface by 
removing the silicon atoms bonded to less than four oxy- 
gens. The oxygen atoms bonded to only one silicon {non 
bridging oxygens) were saturated with acidic hydrogens. 
We note that since the surface of the cavity is rough its 
actual volume V p is unknown and only roughly approxi- 
mated to a lower value by the volume V c of a cylinder of 
radius R c = 20 A. 

Inside the pore a fixed number Nyy of water molecules 
is inserted and the MD simulation is performed at con- 
stant density p ^~ Nw/V c . We assume, according to 
the experimentsaO, that at full hydration the density 
is p f = 0.0297 A" 3 , corresponding to N w = 2600. In 
the following we will present results for N<w = 2600 and 
Nw = 1500 which is roughly the half hydration case. 

Water-water and water-substrate interactions .are re- 
spectively described by the SPC/E site modelo and 
by an empirical potential model, consisting of Lennard- 
Jones and Coulomb forces. The potential parameters are 
those used and tabulated in Ref [l6] where further techni- 
cal details can be found. During the MD runs the atoms 
of the Vycor glass are kept fixed. Periodic boundary 
conditions are applied along the axis of the cylinder (z- 
direction). The motion is confined in the x-y plane. 

A snapshot from the simulation of 1500 water 
molecules is presented in Fig. [jL Only water molecules 
are displayed, the hydrophilic nature of the Vycor surface 
is evident, since all water molecules are attracted toward 
the substrate. 



FIG. 1. Projection on the x-y plane of a snapshot from 
simulation of 1500 molecules of pure SPC/E water inside a 
cylindrical cavity of Vycor glass. The radius of the cavity is 
20 A. For the sake of clarity only water molecules are shown. 
Oxygens are black and hydrogens are gray. 



III. EXCLUDED VOLUME EFFECTS ON THE 
PAIR CORRELATION FUNCTIONS 

The calculation of the site-site radial distribution func- 
tions (SSRDF) in computer simulation is not straight- 
forward when dealing with a confined system, since ex- 
cluded-ii VQ i hw 3£ effects must be carefully taken into ac- 



count .HOE 

From the MD simulation we can calculate in the usual 
way the average number n a L (r) of sites of type j3 lying in 
a spherical shell Av(r) at distance r from a site of type 
a. Then normalizing to the average number of atoms of 
an ideal gas at the same density in the same spherical 
shell, one obtainscj 
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where Np is the total number of (3 sites and V p is the 
volume of the simulation cell. 

At variance with a bulk liquid, if we consider a collec- 
tion of non-interacting particles in a confining volume, 
where periodic boundary conditions are, for instance, 
only along the z-axis, as in our case, the uniform ra- 
dial distribution function is not simply g u (r) = 1. It 
indeed depends on the geometry of the confining system 
and generally will not be a constant. Hence the func- 
tions defined in Eq. (TO) must be further normalized-ip 
the uniform profile g u (r). In our case this is given byt3: 
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Pcyi(Q) is the form factor ofr-the cylindrical simulation 
cell of height L and radius i?j23 
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where j n (%) are the Bessel's functions of order n. 

The properly normalized SSRDF are then obtained as: 
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The correction factor f c accounts for the roughness of the 
surface and the already discussed uncertainty in the vol- 
ume determination. It is adjusted to give the corrected 
pair correlation functions oscillating around 1 at large r. 
The corrected SSRDF for the full hydration case are 
shown in Fig. |2| and are compared to the correspond- 
ing functions of SPC/E water_at ambient conditions. As 
discussed in a previous workED the modifications of the 
oxygen-oxygen (00) and of the oxygen-hydrogen (OH) 
functions relative to bulk water are in good qualitative 
agreement with the experimentsJj For instance the first 
minimum of the 00 function becomes shallower and fills 
in and the H-bond peak becomes less intense upon con- 
finement. 




FIG. 2. Site-site radial distribution functions goo(r) (on 
the left) and goH(r) (on the right) at full hydration at 
T — 298-7^ (full line) compared with the same functions of 
the bulk water (long dashed line) . The functions of the con- 
fined system have been calculated by using Eqs. (PJ)-M). 



IV. LAYERING EFFECT AND STRUCTURE OF 
CONFINED WATER 



Looking at the density profile reported in the inset of 
Fig. H as a function of R 



\/x^ 



y 



2 we see that the 



density profile is very flat in the range < R < 15 A 
and its value is in agreement with the density of confined 
water at full hydration observed in the experiments. For 
R > 15 A water adsorbs and a double layer structure is 
formed close to the surface. The double layer with den- 
sity higher than the average extends for almost 3 A. A 



depletion layer of 2 A, due to the short range repulsion 
of the substrate is also visible. 

We notice that the density profile is almost indepen- 
dent of the temperature. 

As a consequence of the hydrophilic interaction a 
strong distortion of the hydrogen bond network of. w.a-. 
tcr in the layers close to the substrate is observedES'EZI 
Moreover this layering effect has important consequences 
on the dynamical behaviour of water as recently shown.Ej 
It seems therefore appropriate to develop a layer analy- 
sis for the SSRDF by separating the contribution of the 
water molecules which spend most of the time in the in- 
ternal layer (0 < R < 15 A) from the contribution of the 
molecules belonging to the double layer or to the deple- 
tion layer. 



A. Inner layer. 

The SSRDF in the inner layer can be obtained from 
the MD data by using Eqs.( 0)-(|J) where now the radius 
of the cylinder is R c = 15 A. 

In Fig. we show the resulting gooir) for the inner 
layer at different temperatures and compare with the 
equivalent function of ambient temperature SPC/E wa- 
ter. 




FIG. 3. Site-site radial distribution functions goo{r) at full 
hydration calculated for the inner layer (0 < R < 15 A, see 
text) at temperatures T = 298^ (continous line), T = 270K 
(long dashed line) and T = 2AQK (short dashed) compared 
with the corresponding function of the bulk at ambient con- 
ditions (dotted line). The SSRDF of bulk water at ambient 
conditions is almost coincident with that of confined water at 
270K. In the inset the density profile at full hydration for 
T = 298A', T = 270is: and T = 240^ are reported by using 
the same symbols. 

We notice that the goo ( r ) 0I the confined system are 
very similar to those of bulk water. In particular present 
data at T = 270 K superimpose to the oxygen-oxygen 
RDF calculated for bulk water at ambient conditions. 
For the other two SSRDF (not shown) we find the same 
result, suggesting that as far as the microscopic struc- 



ture is concerned water confined in the middle of the 
pore behaves as bulk water at an higher temperature 
(AT ~ 30!^). This result may explain why confined 
water can be more easily supercooled than the bulk. 



B. Outer layer. 

In the calculation of the SSRDF for the outer layer 
the geometry of the confining volume and the form fac- 
tor appearing in Eq.([2j) are different from those used in 
the previous case. Now the function g u (r) must be cal- 
culated in the region 15 < R < 20 A where the water is 
adsorbed on the surface. In doing so, we can assume that 
the uniform density of non interacting atoms is different 
from zero and constant in the range 15 < R < 18 A, 
neglecting the small contribution which comes from par- 
ticles in the depletion layer 18 < R < 20 A. The form 
factor to be used is now that of two concentric cylinders 
of radius R\ — 15 A and R 2 = 18 A respectively: 
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where Vi is the volume of the cylinder of radius Ri and 



PJg(Q) is given by Eq.( |) with R c 



Ri 




FIG. 4. Site-site radial distribution functions goofr) at full 
hydration calculated for the outer layer (15 < R < 20 A) at 
temperatures T = 298K. Thinner solid line is goo(r) defined 
in Eq. (Fy); Dashed line is the uniform profile calculated ac- 
cording to Eq. (El) . The corrected RDF is reported as a thick 
solid line. 

In Fig. |j we report the RDF of the oxygen sites ob- 
tained with this procedure together with the initial dis- 
tribution goo ( r ) calculated according to Eq. (fy) and the 
uniform profile g u {r) for the concentric cylinders. We no- 
tice that the correction for the excluded volume effects 
may still be not completely satisfactory in the region of 
the main peak, due to neglecting the presence of the de- 
pletion layer. In spite of this we can infer that the struc- 
ture of water in the double layer close to the substrate is 
strongly distorted with respect to the bulk. The gooir) 



function does not show any feature typical of the tetrahe- 
dral arrangement of H-bonded molecules. In particular 
the second peak is shifted to higher distances and its in- 
tensity is much reduced with respect to the bulk water 
case. We observe also that the goo( r ) of the outer layer 
does not change much with the temperature, as can be 
seen in Fig. H: only a sharpening of the main peak at the 
lowest temperature is visible. 

The goH{f) and gHH{r) functions calculated for the 
outer layer are reported in Fig. @ at different tempera- 
tures. At variance with the goo\H these functions show 
the main peaks at about the same r-values as their anal- 
ogous functions obtained for bulk water. Yet differences 
are observed in the relative intensity of the first two 
peaks: these may be partially ascribed to residual ex- 
cluded volume effects due to the presence of the depiction 
layer. Also the goo( r ) an d gonir) functions are weakly 
temperature dependent. 




FIG. 5. Site-site radial distribution functions goo{r) at full 
hydration calculated for the outer layer (15<-R<20A)at 
temperatures T = 298 A" (full line), T = 270 K (long dashed 
line), T = 24QK (short dashed line). 




FIG. 6. Site-site radial distribution functions goH(r) (on 
the left) and gnH(r) (on the right) at full hydration calcu- 
lated for the outer layer (15 < R < 20 A) at temperatures 
T = 298K (full line), T = 270.K" (long dashed line), T = 240K 
(short dashed line). 



C. Half hydration. 

The density profile calculated at half hydration (Nw = 
1500) evidence the presence of a double layer of adsorbed 
water molecules, with the same characteristics as in the 
full hydration case. Thus we can calculate the SSRDF 
for the outer layer by using the uniform profile calculated 
in the previous section for the full hydration case (Eq.||) . 
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FIG. 7. Site-site radial distribution functions goo(r) cal- 
culated for the outer layer (15 < R < 20 A) at half hydration 
(full line) compared with the same function at full hydration 
(long dashed line) at temperatures T = 298-K\ In the inset 
the density profile at half hydration for T — 298 A" (full line), 
T = 270A (long dashed line) and T = 240A (short dashed 
line) are reported. 

In Fig. |7j we compare the goo ( r ) functions at the two 
hydration levels investigated, at room temperature. In 
Fig. y the same comparison is done for the goH{r) and 
9hh(t) functions. It is clear from these figures that 
changing the hydration does not affect the SSRDF of 
the outer layer; the same evidence is found at the other 
investigated temperatures. 




FIG. 8. Site-site radial distribution functions gofi(r) (on 
the left) and gnnir) (on the right) at half hydration (full 
line) calculated for the outer layer (15 < R < 20 A) com- 
pared with the same functions at full hydration (long dashed 
line) at temperature T — 298K (full line). 



The SSRDF of the inner layer at half hydration are not 
presented, since the finite volume corrections cannot be 
calculated with enough accuracy when the density drops 
too rapidly (see the inset of Fig. 0) 



V. SUMMARY AND CONCLUSIONS 



We have presented the results of a layer analysis of the 
pair correlation functions of confined water obtained by 
computer simulation. Water is strongly adsorbed, due to 
the hydrophilic interaction with the substrate. Its den- 
sity profile, which is almost independent of the temper- 
ature, shows the presence of a double layer of approxi- 
mately 5 A close to the substrate that we call the outer 
layer. At full hydration water fills the pore and we can 
calculate the structure of the set of molecules in the inner 
layer by excluding the molecules in the outer layer. The 
SSRDF of the molecules in the inner layer are very similar 
to those of bulk water. In particular the structure of the 
confined water in the inner layer at T = 27QK almost 
coincides with the structure of bulk water at ambient 
temperature, thus explaining why confined water can be 
supercooled down to lower temperatures than bulk water. 

On the contrary the microscopic structure of water be- 
longing to the outer layer differs from that of bulk water. 
In particular the oxygen-oxygen pair correlation function 
does not show any signature of the tetrahedral arrange- 
ment characteristic of bulk—water. This is in agreement 
with experimental findingsu and supported by previous 
computer simulation analysis of the distortion of the hy- 
drogen bond network of confined water .113 Moreover the 
SSRDF calculated for the outer layer are almost insensi- 
tive to temperature and hydration level. The first obser- 
vation explains why the temperature derivatives of the 
RDF look so similar to those of bulk waters. The second 
suggests that the structure of water in the inner layer can 
be safely studied in a real diffraction experiment, looking 
at the difference between the data collected at different 
hydration levels. We stress finally that the analysis of 
the dynamical behavior of the molecules confined in the 
outer layer suggests that they_are in a glassy state al- 
ready at ambient temperature. c3 As a consequence their 
relative orientations must be strongly distorted relative 
to the bulk, in agreement with what found in the present 
analysis. The absence of tetrahedral order in the outer 
layer is also responsible for the lower nucleation temper- 
ature upon lowering the hydration level. 
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